A theoretical study of thermal conductivity in single-walled boron nitride nanotubes 
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We perform a theoretical investigation on the thermal conductivity of single-walled boron nitride 
nanotubes (SWBNT) using the kinetic theory. By fitting to the phonon spectrum of boron nitride 
sheet, we develop an efficient and stable Tersoff- derived interatomic potential which is suitable for 
the study of heat transport in sp2 structures. We work out the selection rules for the three-phonon 
process with the help of the helical quantum numbers (k, n) attributed to the symmetry group 
(line group) of the SWBNT. Our calculation shows that the thermal conductivity Atph diverges with 
length as Kph cx: with exponentially decaying /3(r) oc e~'^^'^'' , which results from the competition 
between boundary scattering and three-phonon scattering for flexure modes. We find that the two 
flexure modes of the SWBNT make dominant contribution to the thermal conductivity, because 
their zero frequency locates at k = ±q where a is the rotational angle of the screw symmetry in 
SWBNT. 

PACS numbers: 65.80.-g, 02.20.-a, 63.22.-m, 61.48.-c 
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I. INTRODUCTION 

Recently, experimentalists have accomplished signifi- 
cant progress with the synthesis of boron nitride (BN) 
nano materials, including the hybridized boron nitride 
and grapheneji the boron nitride thin filmj^ and the 
boron nitride nanoribbonj^ Prior to these breakthroughs, 
the BN nano materials have drawn continuous atten- 
tion over past decades. The optical transitions in the 
single-walled boron nitride nanotubes (SWBNT) were 
measured by the optical absorption spectroscopy* or the 
spatially resolved electron energy loss spectroscopy.^ The 
giant Stark effect in SWBNT was observed by the bias 
dependent scanning tunneling microscopy and scanning 
tunneling spectroscopy.^ The lattice dynamics properties 
have been investigated by first-principles calculation^^^?^ 
or the valence force field modeli^i"— For the ther- 
mal transport in SWBNT, both experimental^'''^^ and 
theoretica l "'^^d'' results have confirmed the importance 
of isotopic doping. The theoretical study of the ther- 
mal transport in SWBNT is mainly in the ballistic re- 
gion where Landauer formula is applied with phonon dis- 
persion obtained from lattice dynamics calculations. In 
these studies, the phonon-phonon scattering is ignored 
in the ballistic region, while the effect of phonon-phonon 
scattering on the thermal conductivity in carbon-based 
nano materials has been investigated by the Boltzmann 
equation approachii^r^ However, it is still a blank field 
for the the thermal conductivity of SWBNT. One of the 
main objectives of present paper is to investigate the ef- 
fect of the phonon-phonon scattering on the thermal con- 
ductivity of SWBNT. 

In this paper, we study the effect of the boundary and 
three-phonon scattering on the thermal conductivity in 
SWBNT. In our investigation, we calculate all required 
physical quantities from the Tersoff plus universal force 
field out-of-plane (Tersoff-t-UFFOOP) potential. The pa- 
rameters of this potential are fitted to the phonon spec- 



trum in sp2 BN sheet. It preserves the efficiency and 
stability of the original Tersoff potential. We find the 
selection rules for three-phonon process by using the he- 
lical quantum numbers (k, n) which are assigned to the 
screw and rotation symmetries in the line group of the 
SWBNT. The helical quantum numbers are conserved in 
the three-phonon scattering process and can distinguish 
the flexure modes in SWBNT from other acoustic modes. 
We find that the thermal conductivity of SWBNT shows 
diverging behavior with length , where the 

exponent P{T) decreases exponentially with the increase 
of temperature T . The thermal conductivity is domi- 
nated by the flexure modes, which can be clearly inter- 
preted with the help of the helical quantum numbers. 

The present paper is organized as follows. In Sec. II, 
we flnd the selection rules for the three-phonon scattering 
process. Sec. Ill is devoted to the interaction potential. 
Sec. IV shows formulas for the phonon life time and ther- 
mal conductivity. In Sec. V, calculation results are pre- 
sented and discussed. The paper ends with a summary 
in Sec. VI. 



II. SYMMETRY SELECTION RULES 

In the scattering process, symmetry selection rules im- 
pose strong constraints on the symmetry properties of 
particles which participate in the process. Resulting 
from the scalar property of the scattering operator in the 
three-phonon process^^ the selection rules reflect directly 
the conservation of quantum numbers corresponding to 
all symmetry operations in the system. First of all, the 
energy conservation leads to: 



• + uj' = uj". 



(1) 



We will only exhibit selection rules for combining pro- 
cess. They can be obtained analogously for splitting pro- 
cess and will be listed at the end of this section. Be- 
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sides the energy conservation, the other selection rules 
are governed by the symmetry group of the system. For 
chiral, armchair, and zigzag SWBNT(ni,n2), the symme- 
try groups are the first, fourth, and eighth class of line 
group. The chiral SWBNT have the screw symme- 
tries and the pure rotational symmetries with generators 
S{a, h) and C^. The SWBNT are rotated around z axis 
for a and translated by h in the z direction after the op- 
eration of S{a,h). z is set to be the rotational axis of 
SWBNT. a and h are related to ni and n2.— The rota- 
tional angle for is 2tt/N, where N = gcd(ni,n2) is 
the greatest common divisor of ni and n2. The armchair 
or zigzag SWBNT also have the screw and pure rota- 
tional operations. Besides, there are additional reflection 
symmetries ah/cTv in armchair/zigzag SWBNT. In the 
irreducible representations of the line group, the quan- 
tum numbers corresponding to the screw and rotation 
operations are a set of helical quantum numbers (k, 
with K e {-Tr/h,iT/h] and n = (-7V/2, Using the 
group theory, one finds that the selection rules in chiral 
SWBNT for three-phonon process is merely to guaran- 
tee the conservation of helical quantum numbers (K,n). 
In armchair or zigzag SWBNT, the symmetry property 
is different for phonons in the center or at the edge of 
the Brillouin zone (BZ), i.e (k, n) = (0,0), (±7r//i, 0), or 
(0, ±N/2) for even TV. Those phonons carry the informa- 
tion of the rcfiection symmetries; thus, they arc of higher 
symmetry than the other phonons. From group theory, 
higher symmetry will lead to more complicate selection 
rules. However, those phonons of higher symmetry are 
unable to carry heat energy in the thermal transport, be- 
cause they either correspond to a rigid movement of the 
system or have zero phonon velocity. It should be noted 
that the other long- wave acoustic modes make important 
contribution, except the mode with (K,n) — (0,0). After 
ignoring those phonons, the selection rules in armchair 
or zigzag SWBNT turn to be the same as that of the 
chiral SWBNT. Consequently, the selection rules in all 
SWBNT are: 

K + K = k" mod — , 
h 

n + 71 — n mod N. (2) 

We emphasize that the physical consequence of these se- 
lection rules is to conserve the quantum numbers (k, n) 
corresponding to the screw and pure rotational symme- 
tries of the SWBNT. Besides the simplicity in their se- 
lection rules, (k, n) have several advantages. To the best 
of our knowledge, present paper is the first work to ap- 
ply the helical quantum numbers which correspond to 
the actual symmetries of the SWBNT in the study of 
phonon scattering process in nanotubes. All existing lit- 
eratures use the linear quantum numbers {k,m)f^^ at- 
tributed to the pure translational and rotational opera- 
tions with generators Tqh and C^. q is the number of 
atom pairs in a big translational unit cell and q = q/N. 
The operation Tqh translates the SWBNT by qh in z di- 
rection. All of those pure translational and rotational op- 




FIG. 1: (Color online) Phonon spectrum in BN sheet calcu- 
lated from Tersoff potential, Tersoff-|-UFFOOP potential are 
compared with experimental results from Ref. \33i The UF- 
FOOP potential enhance merely the two out-of-plane (ZA and 
ZO) vibration modes. 

erations together form a subgroup of the line group. Al- 
though Tqh is a symmetry operation in the SWBNT, the 
Cq is not a symmetry operation of SWBNT. As a result, 
the linear quantum numbers [k, m) are not conserved in 
the three-phonon process. The relations between linear 
and helical quantum numbers are n = m mod TV, and 
K — k + TTij- + where j is an integer to keep k in 
its BZ. These relations lead to an additional constraint 
on the selection rules of (fc, to), which has been originally 
pointed out by Dobardzic et al^, and demonstrated by 
Lindsay et al^ 

Analogously, we can obtain the following selection rules 
for splitting process: 

CJ = CJ + CJ , 

K = k' + k" mod ^, (3) 
n = n' + n" mod iV. 

III. INTERATOMIC POTENTIAL 

The interaction for the sp2 BN sheet or SWBNT can 
be described by valence force field models in the linear 
approximationJii^i^ This linear model can be adopted 
to study the ballistic phonon transport without phonon- 
phonon scattering. To investigate the thermal transport 
beyond the ballistic region, it is necessary to apply a more 
realistic interatomic potential which can include the non- 
linear effect and bond reaction phenomenon. Several sets 
of parametric Tersoff bond order potentials have been 
applied in the molecule dynamics simulation of various 
properties for BN materials.—"— However, these param- 
eters are not fitted to the sp2 BN sheet structure. Espe- 
cially, they are not good in the description of the phonon 
spectrum of BN sheet, which is crucial for the investiga- 
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tion of thermal transport in present work. It is important 
to fit a set of Tersoff parameters to the phonon disper- 
sion of sp2 BN sheet, as there is no experimental values 
for nonlinear properties of BN sheet. The fitting proce- 
dure is realized in following steps. Firstly, the phonon 
dispersion in BN sheet is quite similar to that of the 
graphene, because these two materials have similar hon- 
eycomb structure. So we start with the Tersoff parame- 
ters for carbon^ with parameters A and B in the Tersoff 
potential rescaled by a factor of 0.3. Fig. [T] (dash line, 
blue online) shows that this set of Tersoff parameters give 
a good description of the four in-plane vibrations: longi- 
tudinal acoustic (LA), transverse acoustic (TA), longitu- 
dinal optical (LO), and transverse optical (TO) modes. 
However, it results in much lower frequencies for the two 
out-of-plane vibrations: z acoustic (ZA) and z optical 
(ZO) modes. We find that it is impossible to account 
for the in-plane and out-of-plane vibrations simultane- 
ously through adjusting the parameters in Tersoff poten- 
tial where the bond order function depends on the bond 
length and in-plane bond angles. The in-plane transverse 
vibration is captured by the in-plane bond angle term. 
The in-plane longitudinal vibration and out-of-plane vi- 
bration are coupling together in the bond length term. 
As a result, it is impossible to modify Tersoff parameters 
to independently affect the out-of-plane vibration while 
keep in-plane longitudinal vibration unchanged. Brenner 
et al. solved this problem by generalizing the bond or- 
der function to be dihedral angle dependentj2& which can 
treat the out-of-plane vibrations separately. However, 
this is a numerically expensive generalization. Particu- 
larly for large systems, the numerical simulation using 
Brenner potential becomes much slower than the Tersoff 
potential. Actually, this problem can be settled down 
by a more efficient and straightforward approach. We 
introduce a concise potential to manipulate the out-of- 
plane vibration separately. It should be noted that it 
is an attractive way to describe each degree of freedom 
by a separate term, as this can provide clear physical 
picture for each phonon branch. Among various poten- 
tials, we find that the most stable one is the UFFOOP 
potential21: V = Co(Ci + Ci cos(0) + C3 cos (2^)), where 
is the dihedral angle. The four optimized parameters 
Ci are shown in Table HI Fig. [T] (solid line, red online) 
shows that UFFOOP only enhances the out-of-plane vi- 
bration and does not affect the in-plane vibration at all 
in the planar BN sheet. For tubes, especially those with 
small diameters, the dihedral angle depends on the cur- 
vature of the tube, so the out-of-plane and in-plane vi- 
brations are mixing together by the UFFOOP potential. 
The phonon spectrum calculated from Tersoff -I- UFFOOP 
potential agrees well with the experimental values. It 
should be noted however that the theoretical frequencies 
for LO and TO modes do not coincide with experimen- 
tal data in detail, since these two modes are more sen- 
sitive to long-range interactions compared with acoustic 
branchesji2, The long-range interactions originate from 
the charge polarization of boron and nitride atoms in 



TABLE I: Parameters of the four-body UFFOOP potential: 
V = Co(Ci + C2 cos(<^) -I- Cz cos(2<j!))). The second/third lines 
are for configuration with B/N in the center of other three 
N/B atoms. Co is in eV. Other three parameters are dimen- 
sionless. 



bond type 


Co 


Ci 


C2 


C3 


BNNN 


-0.776387 


1.000000 


1.473839 


0.424507 


NBBB 


0.236178 


1.000000 


-0.373869 


-0.223418 



this ionic system. Especially, the crossing of TO and LO 
branches can not be repeated by the Tersoff-|-UFFOOP 
potential. However, we can safely ignore the effect of 
long-range interaction in the study of phonon thermal 
transport, because the LO and TO modes only contribute 
to thermal conductivity through providing new channels 
for phonon-phonon scattering of other phonons due to 
their low phonon velocity. At this point, we would like to 
point out two distinct features of the Tersoff-|-UFFOOP 
potential. Firstly, it is of high efficiency and stability 
which inherits from the original Tersoff potential for car- 
bon system. We have tested that the SWBNT are still 
stable in molecular dynamics simulation at 2000 K af- 
ter 10^*^ steps. Secondly, this potential can give a good 
description simultaneously for the in-plane and out-of- 
plane vibrations of the sp2 BN system; thus suitable for 
the study of thermal transport in these structures. 

IV. PHONON LIFE TIME AND THERMAL 
CONDUCTIVITY 

The thermal conductivity can be calculated by follow- 
ing formula*^ 

'^Ph = X! <nCph(w)wi„,„,^, (4) 
,n,cr 

where Y is the volume of the system. The thickness 
of the SWBNT is considered to be the inter-layer spac- 
ing (3.35 A) in the hexagonal BN multi-layers4 a is the 
polarization index of phonon dispersion. The phonon 
heat capacity is given by Cph — ksx^e^ /{e^ — Vf with 
X = TujJ / {ksT). „ 2 is phonon velocity in the direction 
of thermal current. Under the single mode relaxation 
time approximation, the relaxation rate for phonon mode 
(k, n, a) due to three-phonon scattering process is;^ 

7^ - E 2|^3P^^^.^(A.)iV(.',."),(5) 

where the coefficient IC3P = x ^ x uj'^uj''^ijj"'^ . 

To simplify the notation, we have introduced a single 
prime to denote that the quantity corresponds to mode 
(K',n',(T'). Similarly, a double prime means that the 
quantity corresponds to the mode {k" ,n" ,a"). The two 
modes {K\n',a') and {k" ,n" ,a") are related to each 
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FIG. 2: Fit specularity parameter p to the thermal conduc- 
tivity from molecular dynamics simulation. 



other through the selection rules in Sec. (II). A is the to- 
tal number of atoms in the system. M is the atomic 
mass of a single atom. For the combining process, 
Aw = w + w' - w" and N{uj',lj") = iV,^ - N^' . The dis- 
tribution function is A^o — ~ !)■ For the splitting 
process, Aw = w-w'-w" and A^(w',w") = -I- A^'o + 1 • 
7 is the Griineisen parameter of mode (k, n, a). It is im- 
portant to use mode-specific Gruneisen parameters. This 
has been confirmed by Nika et al. in the calculation of 
thermal conductivity of 2D graphene, which shows that 
the result from mode- dependent Griineisen parameters 
is in better agreement with experiments^ 7 is calculated 
by 7 = —{V/uj){duj/dV) in the 3D structure, and by 
7 — —{S/u}){du}/dS) in 2D system with S as the area. 
In the quasi-one-dimensional SWBNT system, we calcu- 
late it by 7 = — (L/w)(9w/9i). The summation over n' 
in the relaxation rate can be changed into an integral, as 
the SWBNT are usually very long. The relaxation rate 
is then obtained by following formula)^ 



1 



' ps 



where = N^M/ L is the mass per length. Vg — \v' — v"\ 
is the group velocity. The prime over the summation 
indicates the selection rules imposed on those phonon 
modes in the summation. 

In the low temperature region, the thermal conduc- 
tivity is dominated by the boundary and defect scatter- 
ing, as the three-phonon process is very weak. We con- 
sider the boundary scattering by the following relaxation 
ratei^ 



1 

Tbs 



P 



L 



1 



P 



(7) 



where p is the specularity parameter. For example, in 
studies of graphene this specularity parameter was fit- 
ted to experimental data.^*^ Since there is no experiment 



FIG. 3: (Color onhne) Physical quantities in SWBNT(3,3). 
(a). Phonon spectrum parameterized by helical quantum 
numbers (k, n). (b). Phonon spectrum parameterized by lin- 
ear quantum numbers (fc,m). (c). Phonon velocity, (d). 
Griineisen parameter. In (a), (c) and (d), n = -1, 0, 1 are 
displayed by solid line (blue), dashed line (red), and dot line 
(green). The four acoustic phonon modes are highlighted by 
thicker lines. 



on thermal boundary scattering phenomenon, we fit this 
parameter to the molecular dynamics simulation results 
on short SWBNT of 10 nm long where boundary scat- 
tering is important as shown in Fig. [21 The simulation 
is run based on the same potential. In the molecular dy- 
namics simulation, the thermal current is mimicked by 
setting two different temperatures on the two ends of the 
tube. The Nose-Hoove r^'^i^'^ thermostat is employed to 
maintain constant temperatures with a relaxation time 
of 0.4 ps. The description of boundary scattering is phe- 
nomenological, since we do not have an atomic model for 
the boundary scattering. Eq. ^ will give some error in 
the calculation of Vg for short tubes, where the k space is 
discrete. However, the error will be reasonably small, as 
the phonon life time in short tubes is mainly determined 
by the boundary scattering. 

The total phonon life time can be obtained through 
the Matthiessen's rule: 



1 

Ttot 



' ps 



1 

Tbs 



(8) 



From the above, we can calculate the phonon life time 
and thermal conductivity due to boundary scattering and 
three-phonon scattering. This procedure requires the lin- 
ear properties including phonon spectrum w^ „ , phonon 
velocity „, and the nonlinear property of the Griineisen 
parameter 7^„. All of these quantities are obtained from 
the Tersoff-f UFFOOP potential. 
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FIG. 4: (Color online) Thermal conductivity of SWBNT(5,0) 
with length 1, 5, 10 fim from bottom to top. 



CALCULATION RESULTS AND 
DISCUSSION 



We employ the TersofF+UFFOOP potential to calcu- 
late all physical quantities required in the calculation 
of thermal conductivity. Fig. [3] shows the results for 
SWBNT(3,3), where k G (-tt/Zi, tt/Zi] with h = 1.23 A. 
Quantities belonging to n = —1,0, 1 are plotted in solid 
(blue online) , dashed (red online) , and dot (green online) 
lines. Thicker lines are used to mark physical quantities 
of the four acoustic phonon dispersions, i.e., the LA, the 
twisting (TW), and two TA phonon modes (or flexure 
modes). Panel (a) is the phonon spectrum parameter- 
ized by the helical quantum numbers (/t,n). Different 
from the LA and TW modes, the zero point of the two 
TA modes are located at (k, n) = (— a, — 1) and (a, 1). 
Panel (b) shows the phonon spectrum parameterized by 
the linear quantum numbers (fc,m). The BZ of k is only 
half of the BZ for the k and m = -3,-2,-1,0,1,2. 
These curves are obtained through the relationship be- 
tween two sets of quantum numbers. The zero points of 
the four acoustic phonon dispersions are all in the F point 
{k,m) = (0,0). We can not tell the difference between 
the TA modes and the other acoustic modes through 
the position of the zero frequency any more. Panel (c) 
shows the phonon velocity. Let us concentrate on the 
six curves attributed to n = (dashed line, red online). 
Among the six phonon modes at k = 0, LA/TW modes 
are rigid translation/rotation of the SWBNT; so they do 
not contribute to the thermal conductivity. The other 
four modes have zero velocity; thus can not carry heat 
energy during thermal transport. Similarly, the twelve 
phonon modes at k = ±7r//i do not transfer heat energy 
directly, as their velocities are zero. In other SWBNT 
with even N, eg. SWBNT(4,2), we find that phonon ve- 
locities are also zero at BZ boundary {K,n) = {0,N/2). 
Actually, the zero velocity of phonon modes at BZ edge 
and optical modes at F point is the result of the con- 




200 400 600 800 1000 
T(K) 



Q. 



100 



0.01 



10' 
10" 





TA -- 






LA - 








(b) 





200 400 600 800 1000 
T(K) 



FIG. 5: (Color online) (a). Thermal conductivity of 1 
fim SWBNT(3,3) contributed from different acoustic phonon 
modes, (b). Life times for acoustic phonon modes. 



tinuity and inversion symmetry of phonon velocities for 
optical modes. This is also true in bulk materials where 
the selection rules of three-phonon process are merely to 
conserve the momentum besides the energy conservation, 
while the conservation of quantum numbers correspond- 
ing to other point symmetries is ignored. Panel (d) shows 
the Griineisen parameter. Most phonon modes have pos- 
itive 7, except the two TA modes whose 7 are large neg- 
ative numbers. The large negative Griineisen parameter 
was also found in 2D graphene by Mounet et al. doing 
first-principles calculation)^ It is the origin of thermal 
contraction at low temperatures in many materials. The 
obtained Griineisen parameters are quite different for dif- 
ferent phonon modes; so it is important to use a mode- 
dependent Griineisen parameter 7^^ in the calculation 
of thermal conductivity from Eqs. (|4|) ([6]). 

Figure. |4] shows the thermal conductivity in temper- 
ature range [1.0, 1000.0] K for SWBNT(5,0) with sev- 
eral different lengths. In the low temperature region, the 
thermal conductivity is limited mainly by the boundary 
scattering. Kph increases with increasing temperature ac- 
cording to the temperature dependence of phonon heat 
capacity as the boundary scattering is temperature in- 
dependent. Around Tc — 80 K, thermal conductivity 
reaches a maximum value, where boundary scattering 
and three-phonon scattering counterbalance each other. 
Above Tc, the three-phonon process becomes more im- 
portant leading to the decrease of Kph with further tem- 
perature increase. The critical temperature Tc is lower 
in longer tubes where boundary scattering contributes 
less than shorter tubes. With the increase of length, Kph 



FIG. 6: The vibrational morphology of four flexure modes in SWBNT(5,0) with frequency increased from 1.7 cm ^ to 26 cm 
in panels (a) to (d). The values displayed are (k, n,(T, oj) for each mode. 



increases due to weaker boundary scattering, exhibiting 
the ballistic characteristic thermal transport. 

The acoustic modes are major heat carriers, because 
of their higher velocities compared with optical modes. 
Fig. [5] (a) shows the contribution from the four acoustic 
modes to Kph in 1 ^m long SWBNT(3,3). It shows that 
the two TA modes dominant the thermal conductivity, 
although their velocities are smaller than the LA and TW 
acoustic modes. The TA vibrations are also most visible 
in the molecular dynamics simulation. Fig. |6] displays 
the vibrational morphology of four selected TA modes. 
The length of arrow is proportional to the amplitude of 
vibrational displacement. 

The TA modes are more effectively excited than LA 
and TW modes in whole temperature range because of 
their quadratic phonon spectra. As a result, there are 
more TA phonons to carry the heat energy. However, 
the TA modes should not have such a dominant contri- 
bution to Kph if they can not carry heat energy for a 
long time. Hence, it is also necessary to compare the 
life time of different phonon modes. The comparison 
is shown in Fig. [5] (b). It shows that the life time of 
TA mode is about three to four orders longer than the 
LA/TW modes. From Fig. [3] (d), the value of Griineisen 
parameters for general TA modes are close to that of the 
LA/TW modes. It means that there should be almost 
the same phonon-phonon scattering for all three acoustic 
modes between phonon modes (K,n,tT), («;', n', tr'), and 
{k" ,n" ,a"), satisfying the selection rules. However, for 
long-wave LA/TW modes, the energy conservation will 
always lead to a' = a", while the other selection rules 
require k' to be close to k" and n' = n" . From Fig.[3](c), 
two modes with closer k and the same (n, a) will have 
similar phonon velocity, which results in very small value 
of group velocity Vg in Eq. ([B]). As a result, the three- 
phonon scattering rate is extremely strong for LA/TW 
modes. The situation is quite different for TA modes, 
where the selection rules require k' to be far away from 
At" and n ^ n' . It will lead to a large value of Vg , which 
eventually results in a considerably weak three-phonon 
scattering rate. Considering the above two aspects of 
TA modes, we learn that the TA modes are fully excited 
in whole temperature range and can carry heat for a long 
time; so they dominate the thermal conductivity. 

It has been well established that the thermal conduc- 



tivity in bulk materials will converge with length in- 
creasing at some point, yet it is still an open issue in 
low-dimensional materials such as graphene or nanotube. 
The thermal conductivity shows different length depen- 
dence in different transport region. In the pure ballistic 
region, the thermal conductivity is proportional to the 
length; while Kph is length-independent in pure diffusive 
region. Intuited by these two situations, it is a usual trick 
to fit the thermal conductivity as an exponential function 
of length. This technique is useful in the study of ther- 
mal conductivity in nano-materials, where the Fourier 
law fails and the thermal transport is neither ballistic 
nor diffusive. From this point of view, the beta expo- 
nent in length dependence of thermal conductivity is a 
nice way to observe the thermal transport nature. The 
origins for the exponential behavior in ballistic and dif- 
fusive regions are clear, while there is no good reason for 
the exponential behavior of thermal transport between 
ballistic and diffusive. It can be taken as an artifact of 
fitting approximation at this moment. 

We calculate Kph in [1.0, 1000] K for SWBNT(5,0) with 
lengths distributed in [0.3, 10] ^lm and SWBNT(10,0) 
with length L g[0.3, 5] /im. At each temperature, we 
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FIG. 7: (Color onhne) A natural logarithm of the thermal 
conductivity and length. Power function fitting for the ther- 
mal conductivity of SWBNT(5,0) with L G [0.3, lOj^im, and 
SWBNT(10,0) with L G [0.3, 5]^m at 300 K. 
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FIG. 8: (Color online) Temperature dependence for the power 
factors of SWBNT(5,0) and (10,0). 

fit Kph as a power function of L: Kph — aL^ . Two fit- 
ting results at 300 K are illustrated in Fig. [71 In bulk 
materials, = 1 for ballistic transport, while /3 = in 
case of diffusive transport. For nano materials, the ballis- 
tic transport is easier to be observed, because the struc- 
ture is too small for phonon-phonon scattering to take 
effect. However, a pure diffusive thermal transport has 
seldom been seen in quasi-one-dimensional systems42r— 
(3 should be smaller at higher temperatures where the 
three-phonon process becomes more important. A de- 
tailed study of the temperature-dependence of /3 is very 
difficult to be done by molecular dynamics simulation be- 
cause of large computation requirement. The above ap- 
proach in Sec. IV can be used to perform this calculation 
efficiently. Fig. |S] shows that (3 decreases exponentially 
with increasing temperature. j3 has large value at low 
temperatures as it is almost in the ballistic transport re- 
gion. With temperature increase, (3 decreases faster in 
thinner tubes. In the high temperature limit, /? = 0.44 
and 0.42 for SWBNT(5,0) and (10,0). The thermal trans- 
port in both tubes are not in the pure diffusive region. 
Thicker tubes are closer to the pure diffusive transport 
with smaller /3. It is still unclear why the exponent /3 
decays exponentially with temperature increasing, even 
though such behavior has been clearly observed here. 

The divergence of thermal conductivity with nonzero 
j3 is the result of the competition between three-phonon 
scattering and boundary scattering of the TA modes, 
which dominate the thermal conductivity. Fig. [HI shows 
the life time for TA modes at room temperature. The 
drop-off of life time for low frequency phonons is due to 
the quadratic dispersion of the flexure modes, which leads 
to extremely large Griineisen parameter of these modes 
around k = ±a as displayed by Fig. [3] (d). Another im- 
portant feature of Fig. [S] is that the three-phonon scat- 
tering for TA modes around 50 and 70 cm~^ is so weak 
that the life times of these modes are mainly limited by 
the boundary scattering. This part of contribution to the 
thermal conductivity will diverge with increasing length. 



FIG. 9: (Color online) Phonon life time of TA modes in 2^m 
SWBNT(10,0) at 300 K. The life time due to boundary scat- 
tering, Tbs , is shown by blue stars. The three-phonon scatter- 
ing life time, Tpa, is displayed by red open circles. The total 
life time is shownn by black diamond. 

i.e /3 should be 1.0. Other flexure modes have consid- 
erable three-phonon scattering. Contribution from these 
phonons to the thermal conductivity does not depend on 
the length of SWBNT, i.e /3 should be 0. The counter- 
action between these two mechanisms results in a power 
factor p in the range [0, 1]. Actually, a similar phe- 
nomenon was also found by Nika et.al for the thermal 
conductivity in graphene, where the boundary scatter- 
ing still plays an important role for very large piece of 
graphene sample and eventually leads to the increase of 
thermal conductivity with increasing sizej^ 



VI. CONCLUSION 

To conclude, the present work calculates the phonon 
life time due to boundary scattering and three-phonon 
scattering process in SWBNT. The linear and nonlin- 
ear physical quantities required in the calculation are 
obtained from the Tcrsoff-|-UFFOOP inter-atomic poten- 
tial, which inherits the efficiency and stability of the orig- 
inal Tersoff potential and is suitable for the field of heat 
transport. The selection rules for three-phonon process 
are figured out by analyzing the symmetry group (line 
group) of SWBNT. A set of helical quantum numbers 
(k, n) corresponding to the line group is accepted in the 
selection rules instead of the usual linear quantum num- 
bers (fc, m) corresponding to a subgroup of the line group. 
The calculation is focusing on the thermal conductivity 
for SWBNT with different lengths and the contribution 
of different phonon modes in the heat transport. 
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